Controlling electron-phonon interactions in graphene at ultra high carrier densities 
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We report on the temperature dependent electron transport in graphene at different carrier densi- 
ties n. Employing an electrolytic gate, we demonstrate that n can be adjusted up to 4x 10 14 cm -2 for 
both electrons and holes. The measured sample resistivity p increases linearly with temperature T 
in the high temperature limit, indicating that a quasi-classical phonon distribution is responsible for 
the electron scattering. As T decreases, the resistivity decreases more rapidly following piT) ~ T 4 . 
This low temperature behavior can be described by a Bloch-Griineisen model taking into account the 
quantum distribution of the 2-dimensional acoustic phonons in graphene. We map out the density 
dependence of the characteristic temperature Qbg defining the cross-over between the two distinct 
regimes, and show, that for all n, p(T) scales as a universal function of the normalized temperature 
T/Qbg- 

PACS numbers: 73.63.b, 73.22.1, 73.23.b 



At finite temperatures electrons in typical conductors 
are scattered by phonons, producing a finite, temper- 
ature dependent resistivity p [l| . If the temperature 
T is comparable to or larger than the Debye tempera- 
ture 9d — the representative temperature scale for the 
highest phonon energies — all phonon modes are popu- 
lated. In this high temperature regime, p(T) ~ T, reflect- 
ing a classical equipartition distribution of the phonons. 
As T decreases below Qp, however, the bosonic nature 
of the phonons becomes important : only the acoustic 
phonon modes within the phonon sphere of diameter 
k p h = kgT/hv s < kp (where kp is the radius of the 
Debye sphere and v s is the sound velocity) are populated 
appreciably, leading to a more rapid decrease of the resis- 
tivity, p(T) ~ T 5 , known as the Bloch-Griineisen (BG) 
regime [2j-|4j. 

Due to the quasi-elasticity of the electron-phonon (e- 
ph) interactions, the maximal phonon momentum in an 
e-ph scattering event is limited to 2hkp, representing a 
full backscattering of the electrons across the Fermi sur- 
face of radius kp. Since in typical 3-dimensional (3D) 
metals kp is of the size of the Brillouin Zone (BZ), 
2k p > kp, all populated phonons can scatter off elec- 
trons. For low density electron systems, however, the 
Fermi surface can be substantially smaller than the size 
of the BZ, and hence kp -C kp. In this case, only 
a small fraction of the acoustic phonons with energies 
hv s kph < 2hv s kp can scatter off electrons. This phase 
space restriction defines a new characteristic temperature 
scale for the low density e-ph scattering, the BG temper- 
ature Qbg — 2hv s k F /k B < Q D . It was explicitly shown 
in low density 2-dimensional (2D) electron gases formed 
in semiconductor hetero junctions |5J that p(T) drops pre- 
cipitously at temperatures below T < Qbg rather than 
Qp. However, the 3D nature of the phonons in the host 
material, and the low level of tunability of n, make this 
system ineligible for studying 2D BG physics and the 
explicit density dependence of Qbg- 



The advance of graphene jf| 0] brings new aspects to 
the study of e-ph interactions in a low dimensional sys- 
tem. Graphene has a Debye temperature Qd ~ 2300K 
almost an order of magnitude higher than for typical 
metallic systems, and the electrostatic tunability of kp oc 
y/n allows for a wide range of control of Qbg- In addi- 
tion, the single atomic plane structure of graphene pro- 
vides not only a strictly 2D electronic system, but a 2D 
acoustic phonon system as well. These unique properties 
have been considered theoretically, leading to the predic- 
tion that graphene exhibits a smooth crossover behavior 
between the high temperature p(T) ~ T and the low tem- 
perature p(T) ~ T 4 dependence [8(. The slower reduction 
of p(T) at low T as compared to the T 5 dependence ob- 
served in typical 3D conductors can be understood by 
the reduced spatial dimensionality. 

In this letter, we report the experimental observation 
of a 2D BG behavior in graphene. Using an electrolytic 
gate, we achieve extremely high carrier densities up to 
4xl0 14 cm -2 for both electrons and holes, tuning Qbg to 
values of up to -1000 K. In the low T limit, T < Qbg, 
a p{T) — T A is observed, reflecting the 2D nature of 
the electrons and the acoustic phonons in graphene. At 
high temperatures, the resistivity shows a semiclassi- 
cal p(T) — T behavior. From analysis of the experi- 
mental data of p(T), we obtain QBG( n ) an d show that 
p(T) scales as a universal function p(T/Qbg) of the nor- 
malized temperature T/Qbg for all densities |n|, which 
agrees well with the theoretical expressions. 

In previous graphene experiments employing thermally 
grown Si02 layers as the gate dielectric p{T) 
could be measured only in the density range \n\ < 
5x 10 12 cm~ 2 . In this range of carrier densities, p(T) — T 
was reported for all T < 150 K, and, at higher tempera- 
tures, p(T) exhibited a rapid increase presumably due to 
the scattering by thermally activated Si02 polar optical 
phonons [l(| , by thermally quenched graphene ripples [§] , 
or by Coulomb impurities [ll]]. These extrinsic effects 
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FIG. 1. (a) Resistivity as a function of applied electrolyte 
gate voltage V eg at T — 300 K. Right inset: a schematic view 
of the electrolyte gated device. The Debye Layers are formed 
d ~ 1 nm above the graphene surface. Left inset shows an 
optical microscope image (false colors) of a typical etched Hall 
bar device, (b) Inset shows the Hall voltage Vh as a function 
of magnetic field B for different V eg (current across sample 
I =100 nA). The main panel shows the extracted densities n 
from the various Hall measurements as a function of V eg for 
4 different samples. The slope of the line fit represents the 
capacitive coupling of the electrolyte gate to the graphene. 



become less pronounced at higher carrier densities where 
the carrier screening is enhanced [l0| . Furthermore, an 
increased carrier density would result in an increase of 
Qbg-> allowing to access the non-linear p(T) of the BG 
regime in a much wider temperature range. 

In order to increase the electrostatic doping, we employ 
a solid polymer electrolyte gate 12, 13| . This approach 
has previously been applied to graphene samples [14H16| 
and carrier densities of n ~ 10 13 cm" 2 have been re- 
ported under ambient conditions. Following this exper- 
imental approach, we increase the efficiency of the elec- 
trolyte gate by employing a rapid cooling method which 
prevents sample degradation, and reach n > 10 14 cm 
for both electron and holes. 



Fig. 1 (a) shows a working principle of the solid poly- 
mer electrolyte gate used in our experiment. Li + and 
CIO J ions are mobile within a solid "mesh" formed by 
the polymer poly(ethylene)oxide (PEO). By applying a 
voltage V eg to the electrolyte gate, the ions form Debye 
layers on top of the graphene and the gate electrode, re- 
spectively. The extreme proximity of these charged lay- 



ers, separated only by the Debye length A_d ~ 1 nm from 
the graphene surface, results in huge capacitances per 
unit area C eg = 1/eXu- Under ambient conditions the 
electrolyte gate V eg can be swept continuously, doping 
the graphene samples to either electrons or holes, induc- 
ing a modulation of p(V eg ) Fig. 1(b). 

The maximal \n\ that can be induced by the solid 
polymer electrolyte is mainly limited by the on-set of 
electrochemical reactions of the ions with the graphene, 
which typically turn on when V eg > 3 V are applied. 
Although the threshold of the electrochemistry, signaled 
by a steady increase of p with time, depends on the de- 
tails of the particular device and sample quality, the time 
span until complete degradation of the sample is typically 
just a few minutes. We could apply extremely high V eg 
of up to 15 V to the electrolyte, avoiding electrochemi- 
cally induced sample degradation by immediate cooling 
of the sample (< 1 min) below T < 250 K. At this tem- 
perature, both, the Li + and CIO4 ions "freeze" out and 
are no longer mobile within the PEO, fixing the accumu- 
lated charges in the Debye layers to the graphene surface. 
The induced charge carrier densities do not vary signifi- 
cantly over time and temperature until after the sample 
is warmed up again. Accumulated \n\ for each applied 
gate voltage V eg are characterized directly by performing 
Hall measurements (Fig. 1(b) inset). Fig.l (b) shows the 
measured n as a function of V eg . The capacitive coupling 
of the electrolyte can be then estimated from the slope 
of n(V eg ) where we obtain C eg ~ 3.2 /iFcm -2 , which is 
more than 250 times higher than that for typical 300 nm 
thick Si02/Si back gates. 

In this experiment we have measured p(T) for more 
than 10 single layer graphene samples, in the tempera- 
ture range 1.5 < T < 300 K and for \n\ < 2 x 10 14 cm" 2 . 
Fig. 2 (a) shows a representative data set for the mea- 
sured p(T) at various fixed n. Generally, p(T) decreases 
monotonically as T decreases, saturating to po in the 
low temperature limit 17]. This residual resistance po 
stems from almost temperature independent scattering 
mechanisms, such as static impurity and point defect 
scattering, as discussed in previous studies As 
shown in Fig. 2(b) inset, the corresponding mobility fol- 
lows Pq 1 = (poen) m a + bn [li} with the fitting pa- 
rameters a = 3.3 x 10~ 4 Vs/cm 2 and b = 3.9 x 10" 18 Vs, 
representing long and short range impurity scattering, 
respectively. 

At a first glance, the temperature dependent p(T) can 
be subdivided into two different temperature regimes: (i) 
the high temperature linear T re gim e; and (ii) the low 
temperature non-linear T regime [1 91 ] - This transitional 
trend of p(T) at low temperatures can be better scruti- 
nized by subtracting off p from p(T). Fig. 2(b) displays 
Ap(T) = p(T) — po as a function of T in the logarithmic 
scale. At a given density n, each curve of Ap shows a 
clear transition from a linear high temperature behav- 
ior {p r~j T) to a superlinear (p ~ T 4 ) behavior at low 
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FIG. 2. (a) Temperature dependence of the resistivity for dif- 
ferent charge carrier densities of sample G8A4.(b) The tem- 
perature dependent part of the resistivity Ap(T) scales as T 4 
in the low T range and smoothly crosses-over into a linear T 
dependence at higher T. Dashed lines represent fits to the 
linear T and T 4 dependency, respectively. Inset shows the 
mobility fj,o at T = 2 K as a function of the density n. Grey 
line is the theoretically expected mobility due to short and 
long range impurity scattering. 



temperatures, as is expected from the BG model applied 
to electron-acoustic phonon scattering in graphene Q. 
The cross-over temperature between these two different 
regimes appears to be higher for higher carrier densities, 
in good accordance with the BG description presented 
above, where Qbg oc \fn . 

We now quantitatively analyze our data in terms of 
the BG model. Considering the e-ph interaction as the 
major source of scattering, the temperature dependent 
resistivity of graphene can be obtained using the Boltz- 
mann transport theory jsj]: 



A P (T) = l D ' AkF 2 MQbg/T), 



(1) 



where Da and p m are the acoustic deformation poten- 
tial and the mass density of graphene, respectively, Vf 
the Fermi velocity and the generalized BG function for 
graphene is given by the integration form: f s {z) = 
Jq zx (eJcZitr —dx. We remark that Eq.Q]is different from 
a typical BG formula for a 3D metal in three points. 
First, the integrand contains x 4 instead of x 5 , reflect- 
ing the 2D nature of electrons and acoustic phonons 
in graphene. Second, the relevant normalized temper- 
ature scale is Qbg instead of Qd, considering the fact 
Qbg < Od in our experimental range. Third, the ab- 
sence of backscattering for the carriers manifests itself 
in the factor yl — x 2 in the integrand, representing the 
chiral nature of the carriers in graphene |20|. 

Taking the two opposite limits of the temperature 
ranges, Eq. [T] can be further approximated to Ap w 7T 
for T > Qbg and Ap w aT 4 for T < Q BG , where the 
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FIG. 3. (a) Scaling of the prefactors a(n) (Ap ~ a(n)T 4 ). 
Data points were obtained from T 4 fits of the p(T) traces 
at different carrier densities and for different samples. The 
dashed line represents a theoretically predicted fit oc \n\ . 
(b) Qbg at different carrier densities (symbols are defined as 
in (a)). The grey line is a fit to the theoretically predicted 
Qbg = 2hv ay /¥nlk B - 



temperature independent proportionality coefficients are 
explicitly given by [8|: 



TiD 2 A k B 



Ae 2 hp 



m tJ s u F 



and 



12C(4)^4 
e 2 h i p rn vlv 2 F 



(nn)' 



(2) 



(3) 



where £ is the Ricmann-Zcta function. 

Here we particularly note that a oc |n| -3 / 2 , while, 7 
is density-independent. Using these properties, we ob- 
tain 7 and a(n) from the experimentally observed Ap 
at fixed n. First, 7 w (0.14 ± 0.01) O/K is estimated 
from the converging high temperature limit (dotted line 
in Fig. 2(b) for example). This value is in reasonable 
agreement with the previous studies 0, [l(| ■ We then es- 
timate a from each p(T) curve at different densities by 
fitting to Ap ~ T 4 . Fig. 3(a) shows the resulting a ver- 
sus I n I in a wide range of experimentally accessible \n\ 
for 4 different samples. A clear trend of a(n) ~ |n|~ 3 / 2 
can be seen (dashed trace in accordance with Eq. \3§ . 

The combination of the two coefficients a(n) and 7 
allows us to compute the ratios of D 2 A jv 2 and D 2 A / v 5 s , re- 
spectively, and thus evaluate the values of Da and v s 
separately. Employing p m — 7.6 x 10~ 7 kg/m 2 , and 
vf = 10 6 m/sec, we find that the average values for 
each parameter are v s = (2.6 ± 0.4) x 10 4 m/sec and 
Da — (25 ± 5) eV, in a good agreement with the values 
from the literature [1, [ljj, [U Q • Using these we can 
now fit each experimental curve Ap(T) by Eq. 1, using 
Qbg as a single fitting parameter. Fig. 3(b) displays the 
experimentally determined Qbg from this fits as a func- 
tion of n for both electrons and holes for all measured 
samples. The obtained Qbg exhibit the predicted \f\n\ 
dependence, explicitly demonstrating our capability to 
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FIG. 4. Universal scaling of the normalized resistivity 
Ap(T) / Ap(£Q bg) as a function of the normalized temper- 
ature T/Obg, explicitly using the constant £ = 0.2. Data 
points correspond to Ap(T) of sample G8A4 at different n and 
Qbg and are normalized with respect to Ap(£Obg) — the re- 
sistivity at T = £6bg- The dashed line trace represents the 
theoretically predicted scaling of / s (6sG/r)//s(C x ) without 
use of fitting parameters. 



tune the BG temperature of up to ^1200 K with the 
solid polymer electrolyte gate. 

We finally discuss the universal scaling of Ap(T). The 
experimental estimation of Qbg allows us now a direct 
test of the scaling behavior of Ap(T) following Eq. [TJ in 
all temperature ranges. Taking T/Qbg as a dimension- 
less parameter, Eq.[T]can be rewritten in a dimensionless 
scaling form: 



Ap(T) _ UQbg/T) 

ap^qbg) m- 1 ) 



where £ is an arbitrary constant setting the normal- 
ization temperature relative to Qbg- Fig- 4 displays the 
universal scaling behavior of Eq.@] converted from Fig. 2. 
Here we simply choose £ = 0.2 to ensure Ap(^Qbg) is 
within the experimentally accessible range. Remarkably, 
each normalized curve of Ap with different n (thus differ- 
ent Qbg) falls on top of the theoretical curve, indicating 
the BG model of e-ph scattering fully explains p(T) not 
only in the low and high temperature limits, but for all 
temperatures. 

In conclusion, using the electrolyte we have achieved 
extremely high carrier densities of up to \n\ = 4 x 
10 14 cm~ 2 in graphene samples. This advancement al- 
lowed us to observe a strictly 2D Bloch-Griineisen be- 
havior in the measured resistivity, exhibiting the linear 
T to superlinear T cross-over, defined by the gate tun- 
able characteristic temperature Qbg- Our quantitative 
analysis of the temperature dependent resistivity shows 



an universal scaling behavior of the normalized resistiv- 
ity p(T) with the normalized temperature T/Qbgi repre- 
senting the 2D nature of the electrons and phonons along 
with the chiral nature of the carriers in graphene. 
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